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PREFACE 


The U.S. Army Coastal Engineering Research Center (CERC) and the U.S. Army 
Engineer Waterways Experiment Station (WES) entered into an exchange research 
agreement with the Danish Hydraulic Institute (DHI), Horsholm, Denmark, in 
March 1980. In this exchange of research the Corps of Engineers received the 
System 21 Mark 8, which is a proprietary research model used for computing 
water elevations and vertically averaged currents, for use at CERC and WES. 

This report discusses use of the model in estimating short wave transformations. 
The work was carried out under CERC's Wave Estimation for Design work unit, 
Coastal Flooding and Storm Protection Research Program, Coastal Engineering 
Area of Civil Works Research and Development. 


The report was prepared by Dr. Jon M. Hubertz, Oceanographer, under the 
supervision of Dr. C.L. Vincent, Chief, Coastal Oceanography Branch, and 
Mr. R.P. Savage, Chief, Research Division. 


Technical Director of CERC was Dr. Robert W. Whalin, P.E., upon publica- 
tion of the report. 


Comments on this publication are invited. 


Approved for publication in accordance with Public Law 166, 79th Congress, 
approved 31 July 1945, as supplemented by Public Law 172, 88th Congress, 
approved 7 November 1963. 


Colonel, Corps of Engineers 
Commander and Director 
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CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) UNITS OF MEASUREMENT 


U.S. customary units of measurement used in this report can be converted to 
metric (SI) units as follows: 


SS LSS ET I ET ET I I IL I EE I I I EEL IT 


Multiply by To obtain 
inches 25.4 millimeters 
2.54 centimeters 
square inches 6.452 square centimeters 
cubic inches 16.39 cubic centimeters 
feet 30.48 centimeters 
0.3048 meters 
square feet 0.0929 Square meters 
cubic feet 0.0283 cubic meters 
yards 0.9144 meters 
Square yards 0.836 square meters 
cubic yards 0.7646 cubic meters 
miles 1.6093 kilometers 
square miles 259.0 hectares 
knots 1.852 kilometers per hour 
acres 0.4047 hectares 
foot-pounds 1.3558 newton meters 
maida bare 1.0197 x 1072 kilograms per square centimeter 
ounces 28.35 grams 
pounds 453.6 grams 
0.4536 kilograms 
ton, long 1.0160 metric tons 
ton, short 0.9072 metric tons 
degrees (angle) 0.01745 radians 
Fahrenheit degrees 5/9 Celsius degrees or Kelvins! 


lt) obtain Celsius (C) temperature readings from Fahrenheit (F) readings, 
use formula: C = (5/9) (F -32). 
To obtain Kelvin (K) readings, use formula: K = (5/9) (F -32) + 273.15. 
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PREDICTION OF NEARSHORE WAVE TRANSFORMATION 


by 
Jon M. Hubertz 


I. INTRODUCTION 


The U.S. Army Coastal Engineering Research Center (CERC) and the U.S. Army 
Engineer Waterways Experiment Station (WES) entered into an exchange research 
agreement with the Danish Hydaulic Institute (DHI), Horsholm, Denmark, in 
March 1980. One of the first benefits to the Corps of Engineers was the trans- 
fer of the System 21 Mark 8 model to CERC and WES. The System 21 Mark 8 is a 
proprietary research model used for computing water elevations and vertically 
averaged currents. This model has been briefly described in the CERCular of 
October 1980 and in U.S. Army, Corps of Engineers, Coastal Engineering Research 
Center (198la, 1981b)!>»?, 


The System 21 Mark 8 model can be used in two modes--the prediction of 
coastal and estuarine long waves, such as tides or storm surges; and the pre- 
diction of short waves, such as waves with periods of about 10 seconds. The 
steps followed in applying the model are generally the same for both types of 
application. 


This report discusses procedures for using the System 21 Mark 8 model to 
predict nearshore short wave transformation and to make field personnel more 
familiar with the type problems for which it can be used. The example discussed 
illustrates the combined effects of refraction, shoaling, reflection, and dif- 
fraction. The material discussed is intermediate in detail compared with that 
discussed in the publications referenced above and the User's Guide, which would 
be needed for an actual application. 


IIT. APPLICATION ANALYSIS 


1. System 21 Mark 8 Model Grid Setup. 


When applying the System 21 Mark 8 model the first step is to determine the 
characteristics of the waves and the mean water depth in the region of interest. 
By using these parameters and Figure 1, it can be determined whether the model 
is applicable. Given the wave height, period, and mean water depth, the ratios 
H/Lo and D/Lg (Lo = g/2m T*) can be plotted on Figure 1. The model applies 
if the plotted point lies within the shaded region. 


The second step is to determine the size of the model area. This is gov- 
erned primarily by two factors. The area should be large enough to include 
any proposed structure and have boundaries far enough away from regions of 


ly.s. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, "Computer 
Program: $21 Coastal Water Level Variations and Flows," CETN-I-4, Fort Belvoir, 
Va., Jan. 1981a. 


2uLS). ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, "Computer 
Programs - Storm Surge," CETN-I-16, Fort Belvoir, Va., Dec. 1981b. 
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Figure 1. The region of application (shaded area) of 
the System 21 Mark 8 model for short waves. 


interest so that the boundary effects will not influence results. However, it 
should not be so large that excessive calculations are required. Certain 
estimates can be made to determine an appropriate region size. One estimate 
is related to the finite-difference grid that subdivides the region. This 
grid is a set of equally spaced straight orthogonal lines with spacing (AX) 

of at least 10 gridlines per wavelength in the region of interest. This sets 
the approximate distance between the gridlines and helps in estimating the 
region size that can be accommodated. 


A second estimate approximates the time step (At) of the model, expressed 
as 


NE SOO wo 71S) A (1) 


(2gh) 0.5 


where g is gravity and h the mean depth. This states that the time step 
could be about 10 to 15 times the Courant stability criteria for an explicit 
model. However, it may not be practical to have a time step this large since 
the time step might be larger than the wave period; therefore, the wave may 
not be resolved in time. 


An estimate of the computer time needed for a System 21 Mark 8 model run 
can be made using the estimates of the AX, At, size of the grid region, and 
the length of the simulation. This is expressed as 


: Mana Git) 
Run time = Tes sao ae (seconds) (2) 


where M is the total number of gridlines in the x direction, N _ the total 
number of gridlines in the y direction, and T the total time simulated in 
seconds. The number 3,000 is the approximate number of computational points 
per second the model processes on an IBM 3033; on a CRAY computer, it would 

be approximately 20,000. 


It is evident that the grid region, grid spacing, time step, wave charac- 
teristics, and depth are interrelated and that certain limits are imposed on 
the grid spacing and the time step for proper model performance. A number of 
parameter variations may be required before an acceptable combination is 
determined. If waves of different heights and periods are used in the model, 
each case should be checked for conformance to limitations. 


2. System 21 Mark 8 Model Topography. 


Once the model grid has been established, certain information must be 
supplied at each gridline. First the elevation of each gridline must be 
specified. This process is automated since the large number of gridlines 
precludes it being done manually. 


If elevations on the base map are randomly spaced, they can be located in 
the grid by the x and y gridline that is closest. Values of x, y, and 
z (elevation) can then be processed by the interpolation program described in 
Herchenroder (1981)3 to obtain elevations along each gridline. The resolution 
of topographic features in the region is directly related to how much infor- 
mation is supplied to the interpolation program. If a topographic feature is 
not defined by spot values, it will not be present in the interpolated values. 


An alternative method to specifying the spot value of an elevation is to 
specify the location of contours of an elevation. By choosing the contour 
interval, control is kept over which features are to be resolved in the region. 
The contours can be digitized to obtain a set of x, y, and z values for 
input into the interpolation program. This program, in turn, supplies the 
depth values along each gridline. 


III. SYSTEM 21 MARK 8 MODEL COMPONENTS 
1. Initialization Unit (IU). 


The IU creates the input to the computational unit (CU) and also checks for 
errors in the data. Further, it allows the user to make changes or additions 
to a set of basic model parameters that define the type run to be made. Ex- 
amples of such parameters would be the number of time steps, the size of a time 
step, and any options used to print and store the results. In addition, it 
allows for the adjustment of parameter fields, such as the initial surface 
level, the frictional bed resistance eddy coefficients (which are related to 
the diffusion of momentum), and the porosity coefficients (which are related 
to energy loss due to flow through porous media). Boundary data are also 
specified through the IU. These data designate which sections of the grid 


3HERCHENRODER, B.E., "Products From Two Computer Programs Which Process 
Digital Bathymetric Data,'’ CETA 81-13, U.S. Army, Corps of Engineers, Coastal 
Engineering Research Center, Fort Belvoir, Va., Oct. 1981. 
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boundaries are open and the type boundary condition to be imposed on each. 

The program can be run in an interactive mode from a terminal, allowing errors 
in the input data to be corrected efficiently and changes made to the 
parameters. 


2. Computational Unit (CU). 


The CU is the program that does the actual calculations of the water level 
and x,y components of the water velocity along each gridline at every time 
step. It is submitted as a batch job with the user having no control as opposed 
to the IU, which is interactive. The CU uses the data prepared by the IU and 
produces a simulation through the time period specified. The printout of water 
level and velocity components can be obtained at specified time increments. 
Results also can be stored at specified time increments for use later in plot- 
ting the water levels and currents in one, two, or three dimensions at a given 
time or versus time at specified points. 


3. Example Application. 


The System 21 Mark 8 model is applied to an area around the CERC Field 
Research Facility at Duck, North Carolina. Contours of the bathymetry in the 
region modeled are shown in Figure 2. A 60 by 130 grid with a uniform spacing 
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Figure 2. Contours of bathymetry near the pier (October 1980 survey). 
The rectangular boundary near the pier is the region modeled. 
The contour interval is 1 meter. 
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of 16.4 feet (5 meters) was prepared for the bordered region near the pier. 
The wave transformation is expected to be greater within this region than it 
is outside it where the bathymetry is more uniform. 


A wave with a period of 14 seconds and height of 6.9 feet (2.1 meters) is 
input that is normal to the seaward boundary. The value for H/Lo and D/L, 
are respectively 0.0068 and 0.01 using a mean water depth of 13.1 feet (4.0 
meters). These values are within the applicable region as shown by the (X) in 
Figure 1. A time step of 1.0 second is chosen to resolve the wave in time, 
giving 14 time steps per wavelength. The grid spacing of 16.4 feet (5 meters) 
gives about 20 points per wavelength in space for shallow water. Equation (2) 
indicates that a computer run of 15 wave periods requires about 9 minutes of 
IBM 3033 computer time. The actual model run required 10 minutes and cost 
$26. 


A plot of the wave height from various positions offshore is shown in 
Figure 3. The measured values are from gages located along the pier, while 
the computed values are from the model along a line adjacent to the pier. The 
two wave heights (observed and predicted) are not directly comparable. The 
observed wave heights are significant wave height estimates from the peak 
energy of a spectrum, while the predicted values are from a monochromatic input. 
This difference is minimized in this example since the observed spectrum is very 
narrow (Fig. 4). Amother difference between the modeled and the observed is 
wave direction. The best estimate of the wave direction from pressure and 
current sensors at the pier end is toward 238° (North being 0°), which is 14° 
from being straight down the pier. 


This report indicates that the System 21 Mark 8 model provides a reasonable 
reproduction of the measured field wave conditions for the case of a nearly 
monochromatic wave spectrum. Tests representing a broad wave spectra by a 
single height and period are less successful because the energy is spread over 
waves that propagate at different speeds, which violate the basic assumptions 
of the model. Thus, the model's use appears best when restricted to comparison 
studies where test cases (different structure configurations, etc.) can be run 
with either idealized linear or cnoidal waves, or where the basic refraction- 
diffraction patterns are to be studied. The model cannot simulate wave break- 
ing; care must be taken to insure that breaking-wave conditions do not occur 
in a simulation. 


The simulation discussed in this report covered about 0.2 square kilometer 
and cost $26 in computer time. If the same grid spacing, depths, and wave 
periods are applicable, the costs per simulation are about $130 per square kilo- 
meter per model run or $336 per square mile on an IBM 3033 computer. Ona 
faster computer the costs may be reduced by a factor of up to 6, although the 
machine costs may be more. Often, the largest cost is in setting up the grids. 
In any case, practical constraints generally will limit the economic usefulness 
of the model to a small area and only a few runs. If a study envisions the 
need for many runs it may be more economical and perhaps more reliable to use 
a physical model. 


A physical model should be used if wave breaking and reformation is a major 
concern to the effort. The numerical model appears most useful for quick 
studies of small areas where only a few alternatives are being considered. It 
also may be useful when an initial study (perhaps during the reconnaissance 
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Figure 3. The observed and predicted wave heights. 
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The spectrum from wave observations at the end of 
the pier, 24 November 1980. Significant wave height 
at peak frequency is used as input to the model. 


stage of a study or early in the phase I planning stage) with the model is 
needed to help determine the cause of a problem to aid in directing future 
field, laboratory, or numerical investigations required to obtain a viable 
engineering solution to the problem. 


The principal value of the System 21 Mark 8 model is the combination of 
refraction and diffraction, providing improved wave calculations relative to 
refraction models which are more widely used. There are hydrodynamic processes 
the model cannot handle, but it is more sophisticated than most other short 
wave numerical wave codes available. It is also well documented and should be 
considered a helpful addition to the other numerical and laboratory tools 
available for use in analyzing coastal engineering problems. 


IV. SUMMARY 


This report has discussed the range of application of the short wave DHI 
System 21 Mark 8 model for various wave conditions and depths. The general 
steps used in applying the model are presented, along with an example of the 
wave transformation, which simultaneously included refraction, diffraction, 
shoaling, and reflection. 


Possible uses of the System 21 Mark 8 model in the short wave mode are: 
the study of wave transformation near small harbors, shore-connected or 
offshore breakwaters, inlets, or in shallow areas of complicated bathymetry. 
Use of the model for short waves in areas larger than 1 square mile (2590 
square kilometers) for normal wave periods is not usually cost effective. In 
these cases, a combination of models, where one model treats the large-scale 
area and another treats the small-scale area, would be more appropriate. 
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